FINAL report 




During the two years that this grant was in effect, we completed a 
large number of the tasks that were proposed in the original three year 
effort. The remaining tasks are being accomplished on the contract at the 
University of Alabama in Huntsville. In particular, we designed and 
constructed PCM encoder capable of storing data on-board into the mass-memory 
in the encoder at up to 12 megabits per second (total memory of 8 megabytes). 
This telemetry system was programed for two successful flights. The first 
flight was launched in March 1983 from Pokerflat, Alaska. The second, the 
Greenland One Campaign, flew in December/January of 1985. The remaining 
flight, the Critical Velocity Effort, will be launched in May of 1986. All 
parts of the electronic system functioned perfectly during both previous 
flights and the detectors performed perfectly. However, in the first flight 
in Pokerflat, Alaska, an electron arm did not deploy for reasons as yet 
unknown. The ion arm deployed perfectly and good data was acquired. 
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The Greenland flights were spectacularly successful in that all 
elements of the detector system, along with co- investigator experiments, 
functioned perfectly. Efforts are now under way to get these results into the 
literature and AGU presentations in the Spring of 1986. In the data analysis 
area, we have pursued the data analysis of the Condor Critical Velocity 
Experiment which has resulted in several papers included in the appendix, as 
well as one PhD thesis. The Condor experiment stimulated a few independant, 
theoretical investigations, as well as our own continued analysis. We are 
hopeful that a recent model of these releases will help to determine 
perameters that we will use on the launch of the Crit I experiment during the 
Spring of the year. These results will also be presented in an invited talk 
at the COSPAR in July, 1986. 


The electronic system we designed consists of a Main Processor which 
controls all elements of the encoder. The I/O periferals may be accessed 
either by the processor or under automated control via a DMA sequencer. The 
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DMA sequencer allows for the very fast burst capture that would not be 
possible with a regular processor. Thus, the encoder can be run in either of 
two modes — one being, without the DMA under total processor control. In this 
case, we are highly software Intensive and require a large software effort to 
implement our storage. This is the system that was employed on both the 
Pokerflat and Greenland flights. The DMA system will be employed during the 
Grit I experiment and will allow for very fast burst capture, in this case, a 
rate of three megabits per second into one megabyte of memory. 

The processor controls the acquisition of analog digital data through 
several subsystems. The first subsystem is the analog-to-digital convertor 

system \diich has two different possibilities. The first possibility is a slow 
acquisition system of some 40 kllosamples per second at 12 bits into 32 

separate channels. The second is a fast system capable of up to 500,000 

kllosamples per second at 12 bits divided into 16 channels . 

Particle data is accumulated into 16 bit accumulators based on the 
Intel 82C54 counter timer chip. These accumulators are high speed CMOS capable 
of counter-rates up to 20 megacounts per second and controlled in an 
accumulation scheme by Main Processor timing. 

The high voltage sweeps are controlled either in linear or logarithmic 
DAC which is controlled by rights either from the DMA ROM or the Main 
Processor. The entire microprocessor scheme and encoder was built to allow 
for very fast burst capture of data and to allow intelligent control of events 
on the spacecraft. An example of this is the Critical Velocity Experiment in 
which the Main Processor will initiate shape charge ignition via radio control 
under a controlled condition determined by the locally measured magnetic field. 

Science Results 


The flight of 38001 was disappointing in that the main scientific goal 
was not achelved because the two spacecraft were not in the .proper 
orientation. In this experiment, with Gene Wescott of Alaska as the Principal 
Investigator, we were trying to place two spacecraft on the same magnetic 
field line at the time of a radial shape-charge ignition. This was not 
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achieved because one spacecraft severely deviated from our planned trajectory. 

Some energization of local Ions was achieved on the radial shape-charge 

ignition that did occur from the main Instrumentated payload. However, 

because of the lack of electron data, our science results were very limited. 
In fact, our main result was the analysis of a previous flight, 34010, very 
much like this flight, where we did report black holes created in the radial 

shape-a charge distribution and particles energized in that event. 

The Greenland flight was spectacularly successful, in that, due to the 
simultaneous flight of our vehicle, 29023, and University California, 

Berkely's vehicle, we did observe simultaneous precipitation of ions from a 
very far radial distance, as well as, acceleration of electrons simultaneous 
with the deceleration of protons. These reports are now being worked into a 
JGR paper. The further use of the very fast burst capture element of our 
experiment showed, for the very first time, the rapid percipitation of bursts 
of electrons where we can measure the dispersion in electrons from electrons 
accelerated a few thousand kilometers above the spacecraft at the nominal 
altitude of auroral acceleration mechanisms. This is another significant 
result which we are investigating. 

In collaboration with Mike Kelley, we also are investigating the 

creation of possible kllometric radiation due to the resonant emission of 
electromagnetic waves along the resonance cone of whistler waves simultaneous 

with changes in the electron distribution function. In the Critical Velocity 

Area, our main result has been several papers which have now been published 

and one soon to be published on the Condor project which occurred before the 

beginning of this contract, but for which we proposed analysis. The most 

significant result in that effort has been to point out the discrepancy 

between our more negative "resul t^ ahd~thos^”of Pf oj ect “PoFcupihe f epor^^^ by 

Haerendal (1982) in which he reported a large ionization efficiency. 

Hopefully, the resolution of these descrepancles will occur in the upcoming 

Crit I launch in May of 1986. 
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In March of 1988, a barium injection sounding rocket experiment (The Star of Lima) was 
conducted to investigate Alfve'n’s critical ionization velocity (CIV) hypothesis in space. Includ- 
ed in the instrumented payload was a UCSD particle detection experiment consisting of five re- 
tarding potential analyzers. Despite conditions that appeared to be optimal for the critical velo- 
city effect, the particle data , in agreement with optical observations, indicates that a fractional 
ionization of only "5 xlO“^ was observed, indicating that the conditions required for the effect 
to occur are still not well understood. However many of the required phenomena associated 
with the CIV effect were observed; in particular a supertherm al electron population was formed 
at the expense of ion drift kinetic energy in the presence of intense electrostatic waves near the 
lower hybrid frequency. The amount of ionization produced is plausibly consistent with the ob- 
served electron flux, but could also be accounted for by residual solar UV at the injection point. 
It is shown based on the UCSD data set that one obvious explanation for the low ionization 
efficiency, namely that the ionizing superthermal electrons may rapidly escape along field lines, 
can be ruled out. 


1. INIRODUCTION AND EXPERIMENT DESCRIPTION 


The critical velocity hypothesis was first advanced by Alfven [1954] in connection with his 
theory on the origin of the solar system. It states that when a neutral gas moves across magnet- 
ic field lines with respect to a plasma with a velocity at least that given by 




where m„ is the neutral mass, and 4>% the ionization potential of the neutrals, a discharge like 
process can occur which results in the anomalously rapid ionization of the neutrals. The critical 
ionization velocity (CIV) effect has been reported under a variety of laboratory conditions [for 
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example, Piel et al., 1980; Venkataramani and Mattoo, 1980; Aznds, 1978; Danielsson and Bren- 
ning, 1975; FahUson, 1961]. The primary evidence for the critical velocity effect occuring in 
space is a barium injection sounding rocket experiment in the Porcupine series, reported by 
Hatrtndtl |1982|. A comprehensive review of the critical ionization velocity effect in space is 
given by Newell [1985]. 

The Star of Lima was a conical shaped charge barium injection experiment conducted from 
Punta Lobos, Peru in March of 1983. The geometry of the experiment is illustrated in Figure 
1. The Taurus- Tomahawk rocket had a mother-daughter payload: the explosive charge was 
separated by spring release from the plasma diagnostics; at detonation the nominal separation 
was about 2.1 km. The conical charge had an opening full angle of SO degrees, and it directed 
the barium beam almost directly across magnetic field lines towards the instrumentation. The 
launch was just before dawn; the detonation altitude of 429km was calculated to occur below 
the 3200 Angstrom ionizing uv cutoff, but with sufficient visible sunlight for ground based 
television observations (calculations by Weseott et al. [1986] indicate, however, that the barium 
cloud was exposed to about l/40th normal full solar uv). Papers by Kelley et al., jl986j; and 
Weseott et al., [1986] in this issue discuss the electric and magnetic field data and optical obser- 
vations respectively. The present work is concerned with the results of the UCSD energetic 
particle detection experiment. 

The energetic particle detector consisted of five retarding potential analyzers (RPAs), a 
description of which can be found in Carlson [1974]. Two sweeping detectors were devoted to 
ions, two sweeping detectors were devoted to electrons, and one detector measured fixed 
energy electrons. Table I gives the parameters of interest in the energy sweeps. These detec- 
tors exhibit a narrow conical angular response wiUi a half-angle at half maximum of about 3 
degrees. The time needed to complete a measurement of an energy spectrum — 19 ms for the 
fast electron detector, and 38 ms for ail others — afforded better time resolution than most pre- 
vious similar particle detection packages. 

One of the sweeping electron detectors, and one of the sweeping ion detectors, were tilted at 
an 8.5 degree angle to the spin axis to give some limited pitch angle information. The three 
remaining detectors were all oriented along the spin axis, which was vertical to within 10 
degrees; and hence nearly perpendicular to the magnetic field lines (refer to Figure l). These 
detectors thus had their within 10 degrees of anti-parallel to the ram direction of the barium 
neutral beam. 

Elconomic constraints dictated a simple, standard design with a comparatively modest dynamic 
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range of distinguishable particle fluxes. Had the critical velocity effect produced the predicted 
levels of ionization, all five detectors would have saturated, thereby leaving us with a lower 
bound on the fluxes. The much lower fluxes observed are reported in the next section. 

The plasma parameters of primary interest in the Star of Lima experiment are listed in Table 
II, along with various other associated quantities which will be useful in later calculations. 

n. RESULTS: DATA PRESENTATION 

The Taurus-Tomahawk rocket overperformed, so that the solar grazing altitude was 41.5km, 
resulting in the barium cloud being exposed to l/40th the normal ionizing efficiency of full uv. 
The total amount of vaporized barium released in the burst was of the order 10** atoms; the 
amount of barium ionized is estimated from the brightness of the Ball 4554A line observed by 
ground stations to be 5x10*“ ions ( Weseott et al., 1985). This is about l/40th the amount nor- 
mally observed in a similar release exposed to full sunlight. The uncertainties are large enough 
that the precise agreement must be accidental, but it is clear that solar uv is at least potentially 
capable of accounting for all the ionization that occurred. However, it is also true that electron 
impact ionization, in accordance with the critical velocity effect, could have caused as much 
ionization, as is shown below. Moreover, a wealth of plasma phenomena was observed, giving 
insight into the processes in effect. 

Figure 2 shows the electron data near the explosion, which occurred near 340.65 s flight time. 
The Cornell quasi-dc electric field (single-axis double-probe, not despun) measurements are 
included for reference. The ion data is given in Figure 3, which is on the same time scale, and 
includes pitch angle information. Note that the **90 degree" electron sweepihg'“detector 
saturated (while sweeping through the lowest energy step) and took apjH’oximately one second 
to recover. All other detectors appeared to always have dynamic range in reserve; some, partic- 
ularly the "80 ion rpa" (tilted at an 8.5 degree angle to the spin axis) obviously never came near 
saturation. 'Oius the flux rates reported here are believed to be reliable to within a factor of 
about 3 (the uncertainty arising mostly from uncertainty in the energy acceptance widths). A 
chronology of the effects seen by the m aitu diagnostics is given first to provide an overview. 

The ground based television cameras can determine the explosion time to about a single 
frame (1/60 s) resolution. Unfortunately, the time recording on the tapes was hand calibrated, 
so that it cannot be used for a reliable determination of the explosion time. The highest gain 
channel on the Cornell electric field experiment provides the best estimate of the detonation 
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time. Prior to the arrival of any particles, with their associated strong wave activity, a 10 ms 
precursor electric field pulse of amplitude 50/iV/m was detected. This is very likely a magneto- 
sonic wave, which — moving at the Alfve'n speed — would be the fastest traveling plasma dis- 
turbance created by the explosion. This precursor wave arrived at the payload at flight time 

540.658 8. The time of flight of a magnetosonic wave to the payload would only be *3 ms. It is 
not clear how long after detonation such a wave would be emitted; however we may take 

340.658 as the approximate explosion time. For convenience, events in the Star of Lima exper- 
iment will be hereafter be referred to in burst time, where the origin in burst time corresponds 
to 340.658 flight time (21 March 1983 09;58:24.658 UT). 

Figure 4 shows .6 sec of energetic ion and electron data following the explosion on an 
expanded time scale. The analyzing sweep voltage is included in the top panel, with the lowest 
energies in the electrons in panel 2 corresponding to minima in this waveform. ( Refer to Table 
I for sweep ranges). Starting at '62 ms burst time, the first superthermal electrons were 
observed, along with the onset of the oxygen lower hybrid waves (not shown) as rep>orted by 
Kelley et al. [this issue]. Note that counts appear first only in the lowest end of the energy 
sweep; within a few tens of milliseconds the electron tail population heats up (the flux versus 
energy profile broadens). This is the pattern one would expect for electrons heated by local 
wave instabilities. Electrons energized elsewhere would show the inverse dispersion pattern: 
the most energetic electrons would appear first. At 82 ms burst time the first energetic ions 
were observed. Five milliseconds later, the biggest electric field event of the flight occurred; a 
large dc electric field pulse of at least 400mv/m amplitude. Because the orientation of the elec- 
tric field with respect to the single axis of the electric field boom is unknown, this figure is a 
lower limit on the field strength. ~ “ 

To have reached the payload by this time, neutrals would be required to have a velocity of a 
little greater that 20 km/sec, which, in fact, corresponds to an electric field of 440 mv/meter, 
consistent with that observed. However, Weieott et al., [this issue] report no barium ions with 
velocities greater that 14 km /sec. In fact, using their velocity profile, one can show that the 
barium density begins to rise at the payload at about 160 ms burst time ^d reaches a peak of 
about , 2xl0*/cm* at 220 ms, all after this main electric field pulse, but corresponding precisely 
to the times of maximum barium lower hybrid, as later described. IVesumably,then, this dc 
pulse is associated with the arrival of either ambient neutrals pushed ahead of the main barium 
jet or some contaminent of the barium which has never been seen before. To our knowledge, 
no diagnostic payload has ever been in the main jet of such a shaped charge. It is indicative that 
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the primary wave spectrum at this time is more characteristic of oxygen than barium. In any 
event, these early waves do heat the electrons and begin to ionize whatever neutrals are 
present. Thb may in turn cause a strong electric field to be set up in the beam dbection to hold 
back the unmagnetized ions and subsequent polarization of region. The observed electric field 
could then be explained. However, it seems odd to us that beam neutrals at this time, undoubt- 
edly of much lower density than at later times, have such a large effect in accelerating the 
entire background population up to 20 km /sec. Without a full three-dimensional measurement 
of the electric field, the interpretation of the pulse thus remains unresolved. 

Starting at about 150 ms burst time, broad band electrostatic oscillations with a high fre- 
quency cutoff at the barium lower hybrid frequency are observed. A power spectrum of these 
oscillations and a detailed discussion is found in KtUty et al., [I986j. For reference, the local 
beam neutral density begins rapidly to rise at 160 ms, exceeds the ambient at 170 ms, reaches a 
peak of 2x10® /cm® at 230 ms, and decays to a low level (1x10^) by 350 ms. As will be dis- 
cussed in detail below, these times and those of the oscillations correspond also to the time that 
the electron energy flux is gaining energy at the expense of the ion kinetic energy. Note the 
rapid appearance of the peak in the electron spectrum at 165 ms in panel 2 of figure 4. By '260 
ms burst time the oscillations have largely died away. After about 300 ms the ions have ceased 
to lose energy, and the electrons to gain energy. Thereafter, the electron population cools, 
while the ions flux remains stable. The electron population remains isotropic until about 400 
ms burst time. Thereafter a strong loss cone distribution forms in the electrons. 

The absolute value of the ion flux rates can be used to estimate a percentage ionization. The 
peak flux rate, occurring at 242 ms burst time, was about 3x10*® ions/cm ®-s-str implying ion 
densities on the order of "5x 10®/cm®. Using the above neutral cloud density, the ionization at“ 
the time the cloud passed over the instrumented payload can can be readily calculed to be a lit- 
tle more than 10“®. 

A 10 km/s barium atom (a typical value for this type shaped charge) had a gyroradi of 650m 
for our value of magnetic field strength (.22gauss). Since the separation froin the explosion 
point was '2km, we would presumably observe little of the initial thermal ionization; or indeed, 
any ionization which occurred more than a km away from the instrumented payload. Thus any 
ionization which occurred after the cloud passed over the payload would not be detected. In 
this light, it is ^propriate to report dn,/dt; in effect the number of ions observed divided by 
the time the cloud spent within a gyroradi of the payload. Taking, say, a 14 km/s leading edge, 
and a 3 km/s trailing edge one obtains dn,/d(A$ lO® /cm ®-sec. This rate of ion production gives 



a total ionization efficiency (a few times 10“*) that is consistent with the optical observations 
and is three orders of magnitude below the value reported for Porcupine [Haerendel,1982] . 

Electron heating, resulting in a superthermal electron tail, has undoubtedly been observed. 
The peak flux in electrons was about 5x 10’®/cm*-sec-str, which leads to an estimate of the den- 
sity of superthermal electrons of '500/cm* with energy greater than about 7 ev. Photoelectrons 
are much too cold to account for the energetic electron population observed; only collective 
plasma processes could produce the heating that occurred. In particular, there is a clear associa- 
tion of the energetic electrons with rapid electric field oscillations. 

The cross section for electrons to ionize barium is crudely a step function, with the value 1.0 
to 1.2xl0~**cm* above 7 ev, and 0 below 7 ev [ Vainshtein, et al., 1972]. The amount of ioniza- 
tion the observed number of superthermal electrons would be ex}>ected to create is readily cal- 
culated: 


on, 

-^= n,n„<TU, (2) 

where <r =1 xl0~’*cm*; u, =2x10* cm /sec is a typical (10 ev) electron velocity; n, is the 
superthermal electron tail density; and n„ is the density of neutral barium. One finds 
dn,/dl=SxlO*/cm*-sec. This is reasonably consistent with (but somewhat lower than) the rate 
of ion production actually observed as the cloud passed over the payload. 

Figures 5 and 6 show ion and electron energy flux spectra during the event. Each spectrum 
consists of a single upsweep of either the 'BO degree" e- or the "90" ion RPA. Of course, other 
spectra are also available and have been made from the downsweep and from other detectors as 
well, but to .facilitate easy comparison — especially of the quantitative values obtained by 
integrating these curves Table 111 ~ it is desirable to minimize extraneous variations. The 
upsweep is preferred for generating spectra to the downsweep since a wider energy range is 
covered. The detectors chosen had the best coverage (for example, the fast electron detector 
died and hence is unsuitable for the present purposes). 

The earliest ion spectrum — the first for which count rates were significant — is contained in 
Figure 5, 78 ms after the explosion. Tl»is first spectrum has the appearance of a drift maxwel- 
lian with a peak around the retarding potential of 25 ev (implying about a SO or S5 ev actual ion 
energy peak). The second spectrum contained in this figure, taken at .155 ms, shows that the 
peak has shifted to lower energies. Since there has been further ionization, there are more ions 
around at all energies. However it is clear that individual ions are losing energy; and the shift 
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is far too rapid to be accounted for by dispersion alone: Since this is accompanied by the rise of 
electrostatic fields and heating of electrons, it is clear that a collective plasma process is under- 
way transferring kinetic energy from the ions to the electrons. 

The next panel of Figure 5 shows that by 232 ms burst time, the ion distribution has relaxed, 
and lost much of the beam-like form that characterized it earlier. By this time the electron 
heating is well under way. Notice again the shift of the distribution towards fewer ions at high 
energy and more at lower energies. 

Piel et al. |1980] have argued for the importance in laboratory experiments of an inhomo- 
geneous ionization front. Although it appears that such a front could have passed over the pay- 
load, beginning at about 150 ms burst time and ending around 260 ms burst time, this time his- 
tory can also be acounted for by the simple variation of the driving beam density as discussed 
above. These time variations are also seen in the abrupt drop in ion density observed between 
the spectrum beginning at 232 ms and the spectrum starting at 309 ms (refer to Figure 5). 
Thereafter the ion distribution is very stable on the time scale of interest here. For example, 
Figure 5 shows that there is essentially no change in the ion flux distribution between 309 ms 
burst time and 693 ms. All intervening ion spectra are essentially identical and so are not 
shown here. 

The electron spectra show less consistency than the ion spectra, but the general trends are 
clear. Shortly after 100 ms a superthermal electron population begins to form. Figure 6 shows 
first the electrons at 232 ms burst time, and reaching their peak at 309 ms. Although the ions 
experience a sharp drop off at this time, the electrons do not; Figure 6 shows the electrons are 
still at peak energy flux as late as 386 ms. Thereafter the electron flux does begin to drop, 
however, and as the final panel of Figure 6 illustrates, the electron energy flux is much reduced 
by 463 ms burst time. 

Some sample energy fluxes and densities calculated by integrating these curves are summar- 
ized in Table III. In computing these values, the ion dbtribution was assumed to be extended 
over one ster-radian, and the electrons over 4v ster-radians (in accordance with the apparent 
isotropy of the electrons as discussed below, and the observed strong anisotropy in the ion dis- 
tribution). 

Figure 7 gives some limited insight into the pitch angle dependence of the energetic electron 
population. The plots in this figure are generated from the RPA which was devoted to elec- 
trons and which was tilted at an 8.5 degree angle to the spin axis. Thus it spun through about a 
17 degree range of pitch angles. The detector was sweeping in energy; therefore the data in 
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Figure 7 is a subset, each panel corresponding to a given energy level. Note that for approxi- 
mately 0.4 seconds the population is apparently isotropic (or at least flat near 90 degrees). This 
is the time of peak electric field fluctuations, and also the time of peak particle fluxes. 
Thereafter the population becomes sharply peaked, as all electrons not very near 90 degrees 
quickly escape along field lines. At these later times, the ratio of the population near 90 to that 
at about 70 is roughly S0:1. We will return later to the significance of the initial isotropic 
phase. 


III. Discussion 


The data shows that the phenomena associated with the critical velocity effect -- namely heat- 
ing of the electrons at the expense of ion kinetic energy, and associated broad band electric 
waves around the lower hybrid frequency -- were indeed observed. The process failed to build 
on itself in a rapid non-linear fashion; i.e., it did not reach discharge. This is something of a 
puzzle when contrasted with the results of Porcupine, a similar barium injection experiment in 
which nearly 30% ionization efficiency was reported [Baerendel, 1982). 

A number of possible reasons for the low ionization efficiency ("5 xl0“‘) observed in the 
present experiment can be advanced. To discuss all of these in any detail is beyond the scope 
of the present work (an article discussing many of the possibilities in detail is in preparation). 
However our data provides considerable evidence concerning one apparent difficulty; the r^id 
escape of the superthermal electrons along the magnetic field lines. 

For the critical velocity efiTect to be viable, it is clearly necessary that a newly energized 
superthermal electron have time to have at least one ionizing collisioh during the course' of the 
process. Thus, in analogy with the Townsend condition for beam plasma discharges, Brtnnmg 
[1982] has introduced the requirement; 


Tp> Tg+Tion 



where fp is the time over which CIV efl'ects occur, Tg v& the energization time for a superther- 
mal electron, and r is the time it takes such an electron to have an ionizing collision with a 
neutral. Because of the rapid time scale for the lower hybrid instability, it is probably safe to 
take the energization time, Tg, to be much shorter than the other times of interest. 
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The time it takes an electron to have an ionizing collision is readily calculated: 


.(*') = ■ 




■fts 25 ms 


(4) 


where the cross section [ Vainthtein ti ai, 1972] and velocity used were those appropriate to a 
10 ev electron; and the neutral barium density was taken to be 2xl0*/cm*. The time available 
for the CIV process to grow to discharge is not well defined, but is of the order of the 'tontact 
time” that fast flowing neutrals spend on a given field line. This is hundreds of milliseconds, 
clearly much longer than the necessary 25. Thus the Townsend condition in this sense appears 
satisfied. 

From equation (2), it is clear that this r^„ is also the time scale for the growth of the plasma 
density, which is essentially independent of how rapidly the electrons are heated. Their energy 
is fixed by a resonance condition with the wave. One could argue that, given only 200 ms of 
'Contact time,” the density at the payload could only increase by a factor of a few hundred. 
However, this would not explain the absence of ionization at, say, 6 km from the explosion, 
where the time scales are 3 times as long. At these distances. Porcupine produced its maximum 
ionization |Haerendel,1982] . 

There is also, however, the question of electron containment. A 10 ev superthermal electron 
has a velocity of 2000 km/s, and hence could escape in a fraction of a millisecond from the 
cloud, which has a length parallel to B of only ‘600 meters. Thus the escape time appears to be 
-much -more rapid than the ionization time, which would suppress CIV effects. Now such an 
escape current would tend very rapidly to build up a containing electric field; the question is 
whether the background plasma — which has an electron density much larger than the superth- 
ermal electron. density — would rapidly short out any such polarization field parallel to B. It will 
now be shown that our experimental data shows that the electrons were indeed trapped for 
much longer than the nominal escape time. 

Refer to Figure 7, which illustrates the dependence of flux rates on pitch angle ^ various 
fixed energies. Notice that during the time of peak flux rates, say the .3 sec starting at burst 
time 92 ms, the electron distribution appears to be independent of pitch angle. This is not, of 
course, to say that the flux rates are constant during this interval: indeed there is time depen- 
dence as is discussed in detail in the previous section. However the fluxes appear to have no 
dependence on pitch angle during this time period. 
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At a somewhat longer time a marked loss cone distribution emerges. Indeed, after say 600 
ms burst time, ^the ratio of the flux neair 90 to that at about 75 degrees is (with appreciable vari- 
ation) around 30:1. This, however, happens on a time scale much longer than CIV processes 
could possibly be operative. Thus the pitch angle data indicates that the hot electrons are con- 
tained, at least for a much longer time scale than the fraction of a millisecond nominal escape 
time. 

It may be suggested that the electrons are being energized primarily along the field lines, as 
well as escaping along them, in such a way that the observed isotropy is produced. It is a priori 
unlikely that these effects could balance so closely. Furthermore, there is a compelling argu- 
ment based on energy considerations which shows explicitly that such a coincidence is not 
occurring. 

Refer to Table III. This table shows the ions initially losing kinetic energy, and the electrons 
heating up. Note that the energy density of the ions is at all times much larger than the energy 
density of the electrons. However the electron energy flux is much higher than the ion energy 
flux; indeed, by burst time 386 ms, the electron energy flux is 2.4xl0*^ev/sec-cm*, whereas the 
ion energy flux is only 5.7 xl0*°ev/sec-cm*. This is not surprising, since for a given energy 
density, the much lighter mass of the electrons implies a much higher energy flux. The point is 
that these differences in energy flux could not be maintained if the electrons were escaping 
along field lines. 

To conserve energy, it is clear that the electron energy flux escaping along field lines can only 
be a fraction of the ion energy flux into the region perpendicular to B. Unfortunately it is not 
possible merely from the observations of total net flux to deconvolve the sinks and sources. To 
be more explicit: the ion flux, for exaihpile, is increased by new ionization, and decreased by 
dispersion as well as by energy transfer from the ions to the electrons. Nonetheless, it is an 
observational fact that large numbers of ion are not being created, which allows a very rough 
estimate of the electron trapping time — which will be denoted rj- — as follows. The physical 
quantity which most closely corresponds to the detector count rate b the electron energy flux, 
jgf. Imagine, a box of width w p>erpendicular to the beam direction and arbitrary lengths in the 
other two dimensions. At any instant of time, the electron energy flux out of the sides of this 
box along B b proportional to . If the electrons are trapped, thb energy loss b reduced by 
the ratio of the un-imp>eded electron transit time to the trapping time. Thb b essen- 
tially the length w divided by a typical electron velocity. Now, -the npp>er limit to the energy 
available in this box b given by the difference between the ion energy flux ( ^,) which b 
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flowing in and that which is flowing out the other side. Setting this value as a maximum, we 
arrive at: 


»■/«« 


( 5 ) 


where the derivative of is taken along the beam direction. The maximum value of this 
derivative is certainly given by dropping the entire value of over one ion gyroradii, about the 
value of w at the location of the spacecraft. Using the values in Table III at S86 ms for and 
jEi, and a value of about 0.3 ms for r/„,, we find that the trapping time is greater than about 
25 ms. This estimate is obviously very crude, but is remarkably consistent with the ionization 
time given in equation 4. In feu:t, rj- is really an energy trapping time and should never be 
greater than the ionization time for weakly ionizing discharges. Thus, the effects of electron 
escape appear to be no greater than the damping effect of energy loss due to the ionization pro- 
cess itself. The present arguments make clear that energy conservation requires a trapping time 
orders of magnitude larger than the nominal fraction of a millbecond escape time the electrons 
would have if they were uncontained. 

Thus both the energetics of the observed particle distributions, and the isotropy of the elec- 
tron distribution demonstrate that the superthermal electron population is contained for far 
longer than the nominal escape time. 

IV. Summary and Conclusions 

An energetic (1-50 ev) particle detection experiment was designed and built for inclusion 
on the Star of Lima, the first sounding rocket experiment intended solely to test the critical 
velocity effect in space. The results of the experiment were negative, in that the ionization 
efficiency from superthermal electron impact was quite low (k5x 10~^ fractional efficiency). 
However, a detailed analysis of the particle data demonstrates some important results. As 
required by all critical ionization velocity theories, superthermal electron tail formation was 
observed. In accordance with predictions by prevailing theories, this was accompanied by 
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broadband electrostatic oscillations near the barium lower hybrid frequency and a demonstrable 
energy loss in the ion spectrum. There b thus no doubt that the electrons were energized by a 
collective plasma process at the expense of the ion drift kinetic energy. Thus some of the 
essential features of the CIV mechanism were directly measured. Analysis of this data puts 
constraints on speculations into why the CIV operated at such a low level. 

The most important example of this is that the superthermal electrons were trapped for far 
longer than their nominal escape time. The pitch angle isotropy observed, as evidenced by Fig- 
ure 7 suggests such a containment. Energy conservation arguements based on the measured 
fluxes demonstrate convincingly that this is so. Thus one important possible explanation for 
the observed low efficiency is eliminated. 

The experiment described in this paper was performed under conditions which would ordi- 
narily be expected to maximize CIV effects, and yet only a very low ionization efficiency was 
observed. It is thus established that the counterstreaming of a plasma and a neutral gas in a 
space environment does not necessarily lead to CIV effects. 

In view of the large number of proposed applications for the critical velocity effect in space 
[Newell, 1985] it is of considerable practical importance — as well as being of intrinsic funda- 
mental interest — to resolve the question of when to expect CIV effects in space. Further mag- 
netospheric barium injection experiments remain the most promising means of resolving the 
issue . 
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Figure Captions 


Figure 1. Geometry illustreU^ing the Star of Lima barium injection experiment. The arrows 
pointing out of the instrumented payload show the direction of the "90" and '180" degree RPAs. 


Figure 2. Electron differential energy fluxes near the explosion (t=0). The Cornell electric 
field data (not despun) is shown for reference. 

Figure 3. Ion differential energy fluxes near the explosion. 

Figure 4. The electron data near the explosion on an expanded time scale. Again Cornell elec- 
tric field data is included for comparison (refer to Kelley, et al., this bsue). 

Figure 5. Ion differential energy flux as a function of energy at selected times. Each curve is 
based on the data from one upsweep of the "90 degree" ion detector (which took 38 ms to com- 
plete). The times given are the burst time at which the sweep was started. 

Figure 6. Electron differential energy fliix“ as a function of energy" at selected times. Each 
curve is based on the data from one upsweep of the "80 degree" electron detector (which took 
38 ms to complete). The times given are the burst time at which the sweep was started. 

Figure 7. Electron differential energy fluxes at specified energies (6-11.4 ev) for the "80 
degree" detector. The actual pitch angle is also plotted. Note the period of ^{H-oximate iso- 
tropy during the time of peak fluxes (the lower energy electrons are becoming more numerous 
during this time period) . ’ ' 
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Table Captions 


Table I. Characterizing the energy sweeps of the retarding potential analyzers. "Accum"is the 
accumulation time at a single energy step, and 'Spectrum” is the time needed to complete one 
upsweep or downsweep. 


Table II. Important parameters and typical values for the Star of Lima experiment. 

Table III. Energy densities and fluxes for superthermal electrons and ions at times near the 
burst. Values are obtained by integrating the area under curves such as those shown in figures 
5 and 6. 
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Table 1. RPA Sweep Parameters 


Particles 

Angle 

Accum (msec) 

Spectrum (msec) 

Energies (ev) 

e- 

■90 " 

2.40 

19.2 

1.5-47 

B 

■ 90 " 

4.80 

-* 

6.5 

B 


4.80 

38.4 

1.5-48 

ions 

"90 " 

4.80 

38.4 

.8-47 

ions 

■ 80 " 

4.80 

38.4 

.8-46 



























Table II. Star of Lima Plasma and Particle Parameters 


Altitude 

428km 

Background plasma density 

n = 2 or 3 X 10*/cm* 

Background neutral density 

n, =6.7xlO’/cm* 

Magnetic field 

B == .22 gauss 

Debye length 

Ax) =2cm 

Typical Ion gyroradius 

p, (Ba, 10km/sec)=650m 

Typical e- gyroradius 

p, (lOev) = 52cm 

Ion gyrotrequency 

(Ba) = 15.3 rad /sec 

Electron gyrofrequency 

u)„ =3.87xl0Vad/sec 

• 

Plasma frequency 

Uf, =9x 10‘rad/sec 

Ba Ion Plasma frequency 

Ufi = 1.8xl0Vad/sec 

Lower hybrid frequency 

Wft =6.5x lO’ rad/sec 

Alfven velocity (0) 

Px = 700 km /sec 

Particle/Field pressure 

pa 2x10’* 































Table m. Selected Energ>’ Densities and Fluxes 


Time 

Flux (10'*ev/sec-cm*) 

Density (kev/cm*) 

msec 
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ions 
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ions 

78 

1.4 

.063 
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123 

155 

3.8 

.45 

17 

1151 

232 

1.2 

.21 

9.5 

792 

309 

m 

.041 

18 

182 

386 

m 

.057 

18 

223 

463 

.14 

.054 
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272 

693 

.33 

.042 

3.2 
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ABSTRACT 
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High explosive shaped charge experiments King Crab, Bubble Machine I 
and II, designed to perturb the ambient plasma and magnetic field were 
flown above 460 km on Taurus Tomahawk rockets from Poker Flat in March 1980,' 
1981 and 1^82^_repectively . The last two flights were a mother-daughter 
combination with the instrumentation section remaining attached to the 
rocket. Cornell University furnished a single axis dipole electric field 
detector, and a fixed bias cylindrical Langmuir probe. A three axis 
attitude magnetometer was also flown. The University of California at San 
Diego furnished a curved plate energetic ion and electron electrostatic 
analyzer. 

The first experiment. King Crab, revealed the existence of a barium 

plasma depleted region or black hole of about 5 km diameter centered on the 

1 

burst. All electric power to the instruments failed on the liftoff of the 



1981 Bubble Machine I, but useful optical data were obtained. The third 
flight a year later produced useful data from all but the electron detector 
on the ESA. Ground based optical and telluric field instrumentation recorded 
evidence for Injection Induced waves of about 5 second period. Delay 
times Indicate a slow (175 - 385 km/sec) propagation from the burst point. 
Concerning locally produced electrostatic Instabilities the wave measurements 
Indicate that finite Larmor radius stabilization occurs for the weakly driven 
low (barium) density situation, but that shorter wavelength waves can occur 
early in the expansion event when the barium density is higher. These 
results are in reasonable agreement with a theory described by Sperling and 
Krall (1981). 
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I. 


Introduction 


Taurus Tomahawk rockets containing high explosive shaped charges with 
barium liners were launched from Poker Flat, Alaska in March 1980, 1981, 1982. 
The shaped charges were designed to make waves by perturbing the ambient plasma 
and magnetic field by injecting a thin sheet of barium gas radially perpen- 
dicular to the B field. The 1981 and 1982 experiments were a mother-daugh- 
ter combination, with the mother containing instrumentation for measuring 
fields and plasma parameters about a km away from the detonation. An 
electrical power failure at lift-off on Bubble Machine I eliminated the 
instrumented experiments, but the repeat experiment Bubble Machine II a 
year later worked well. Ground based data were obtained in all experiments. 

The instrumented section of the Bubble Machine II rocket remained attached 
to the Tomahawk motor and consisted of a single axis (3 m tip to tip) 
dipole electric field detector operated from dc to VLF frequency, and a 
fixed bias cylindrical Langmuir probe supplied by Cornell University. 

A fixed energy (240 eV.) electron retarding potential analyzes, (RPA) and 
curved plate energetic ion and electron electrostatic analyzer (ESA) were 
supplied by the University of California at San Diego (UCSD). One of the 
ESA's was devoted to ions; the other to electrons. Both swept from about 
50 eV/q to 10 keV/q with a spectrum obtained in approximately 125 m sec. 

Good data were obtained from all detectors except for the electron channel 
on the ESA. Space was quite limited and precluded inclusion of a full 
attitude determination system. A three axis magnetometer yielded aspect 
with respect to the magnetic field. It should be noted that a single axis 
dipole detector can only receive a single component of any (electric field 
fluctuation which takes place more quickly than one half of the rocket 
spin period. 


Launches were planned for quiet conditions so that the effects of the 
perturbation would be evident in situ and at remote sensors on the ground. 

A telluric current recorder located at Poker Flat measured the voltage 
across a 1.597 km long, magnetic N-S oriented grounded dipole. VLF emissions 
were also recorded from an antenna located on a ridge near Poker Flats. 

The Injections occurred above 460 km altitude, and were observed with 
optical Instrumentation from sites in Alaska and from the NASA Ames Learjet. 

The shaped charge detonation did make waves. Various aspects of the 
experiments are discussed in the following sections. 



II. Data Presentation 


High Explosive Radial Shaped Charges : Neutral Gas Dynamics The high 
explosive shaped charges used in the Bubble Machine I and II were an improved 
version of the first King Crab device (Wescott et a1 , 1981). Figure 1 
shows a schematic drawing of the device. 4.5 Kg of comp. B high explosive 
formed the shaped charge surrounding the 1.04 Kg V>shaped ring liner of 
barium metal. Detonation occurred at the center of the shaped charge. 

The pressure of the explosion on the barium liner causes It to collapse 
and partially vaporize. We have assumed 151 vaporization, but this number 
Is not well known. Michel (1974) concludes that 27% vaporization Is the 
theoretical limit In a 30* conical Ba shaped charged. Michel, (1969) 
measured 11% In a small laboratory shaped charge and estimated 20% for a 
full scale model. Some of the vapor Is accelerated to high velocity radially 
with a small angular distribution about the plane perpendicular to the 
axis. Figure 2 shows the velocity distribution measured from TV signal 
Intensity for Bubble Machine II. Initially, the barium gas has a small 
Initial fraction of Ions created In the explosion process. Estimates of 
the thermal Ionization fraction range from n ■ 10“3 to < 10*^, Haerendel 

.11 982) . Further Ionization 1 s produced by by several ^processes. The 

releases are carried out In full solar u.v. radiation which produces 
Ions with an exponential time constant of about 20 seconds Michel (1969). 

In the case of radial Injection perpendicular to the field, the gas with 
velocity greater than '2.7 km/sec exceeds the Alfven critical velocity 
for self Ionization, and when the density is appropriate some Ions should 
be produced by the Alfven mechanism (Alfven, 1954; and Arrhenius, 

1975; Deehr et al., 1982; and Haerendel 1982). 
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Figure 2. DIFFcrenttal IJa atom numbers vs radial velocity distribution, 

Oubbic Machine 11 derived from IV signal level scans. Fifteen 

pe r ^^bcr was assumed to assign number scale 



There have been four radial shaped charge injection experiments (Wescott 
et a1., 1980, Deehr et al., 1982 have described the first two). Although 
differences have been observed between them, we can describe the typical 
features. The barium metal liner is vaporized and comes out in a thin 
sheet radially. As soon as details can be detected In TV or film camera 
pictures (about 3 seconds after detonation) the neutral barium gas can be 
seen to have finger like wispy jets at the leading edge, and blob-like 
irregularities in the slower portions. The wispy jets have velocities 
ranging from 6-8 km/sec. The annulus of blob-like irregularities has 
material traveling from 3 to 4.5 km/second. The bulk of the gas has velocity 
from 0.94 to 2.6 km/sec and is in an inner annulus also with irregularities. 
Bubble Machine II was observed from a NASA Lear jet flying over British 
Columbia. The view was almost from the side allowing calculations of the 
thickness of the torus of slow barium gas (Figure 3). When the slow ring 
had expanded to a diameter of 45 km, video signal scans across the major 
and minor axis provided a measure of a half angle of 13 degrees for the 
torus. 

Barium Plasma Effects : Xing Crab Revisited 

As the neutral gas expands outward any barium Ions formed are held by 
.the magnetic fiel d l ines and -left behind — In-the absence of -an electric - 
field as in the King Crab experiment, (Wescott et al., 1980) the ions 
appear in a stellate "wagon wheel" configuration. There is a central 

region depleted in ions which appears as a black hole. For reference 

} 

Figure 4 shows the Ba ions left behind in the "King Crab" experiments. 

The diameter of the black hole was about 5 km. 

The King crab experiment was the first radial shaped charge experiment, 

launched on a Taurus Tomahawk rocket 16 March, 1980 (Wescott et al., 1980). 

Detonation occurred at 1202:10UT at an altitude of 571 km. There was 
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Figure 3. 


Side view of Bubble Machine II from geophysical Institute TV on 
NASA Ames Learjet, 21.5 sec. after detonation. The plane of 
the low velocity Oa gas doughnut is 17" with respect to the 
viewer and the distance is 1950 km. The two stars are v Andromeda 
and B.S. 157. 
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almost no convection electric field as the stellate Ion structure was 
observed for 11 minutes without drift motion and distortion. The Chatanika 
Incoherent scatter radar also detected no electric field (J. Vickery, private 
communication, 1980). However background levels of YLF and ELF noise were 
evident at the burst time. The auroral background at 4278 A was also high 
and pulsating at the burst time. Therefore the natural background noise 
masked any detonation produced effects. 

The black hole and Ion rays are the most Interesting feature of the 
King Crab experiment. There were about 18 bright rays emanating outwards 
from a 5 km diameter region depleted In Ions. Figure 4 shows the TV view 
from Ft. Yukon, Alaska looking nearly In the magnetic lenlth. From Ester 
Dome, about 25** away from the magnetic zenith, the black hole Is elliptical 
with elongation In the magnetic field meridian. Understanding this feature 
was one of the goals of the Bubble Machine I and II experiments. Before 
presenting the new data, we first summarize our knowledge of the black 
hole formation based on the King Crab results. 

The observed velocity distribution of the neutral Ba atoms (Figure 2) 
shows a pronounced peak at 1.6 km/sec of 3.4 x 10^3 atoms/km/sec. At 1 
km/sec. there are still 1.8xl023 atoms/km/sec. At these low velocities 
the -Ion gyro-radius Is -49 and 31 m respectively _for these- two _vel.oc1 ties, ___ .. 
so any 1on created by solar u.v. Irradiation should gyrate closely around 
the field line where It began. If nothing acted to change the situation 
there should be a high concentration , of Ba Ions surrounding the burst 
point. Some efficient process must haveiacted to allow all or roost of the 
barium Ions to escape from a black hole of diameter about 5 km. 

If the magnetic field were reduced to a low level, .the Ions as created 
could continue radially outward as the gyroradius Is proportional to 1/B. 

For Instance a 1.6 km/second Ion In a 960 nT field has a gyroradius of 2.5 km. 
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Figure 4. TV view of the stellate Ion structure produced In the March 16, 
1980 King Crab experiment 205 seconds after detonation. The 
inner hub region depleted of ionization is called a black hole. 
An almost comi)lete absence of an ionospheric E field allowed the 
structure to be seen for about 20 minutes. 
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There is not enough kinetic energy in the neutral Ba gas to create a magnetic 
bubble that large. Assuming 15% vaporization of the liner, and the differential 
velocity distribution of Figure 2, the total kinetic energy is 540 kJ. If 
this energy was totally dissipated in annihilating the magnetic field, a 
spherical magnetic bubble of diameter 1012 m could be created. However, 
we know that the neutral Ba expanded outward from the burst point with no 
apparent retardation, so only a fraction of the kinetic energy was lost. 

In Figure 5 we have plotted the collision frequency between a magnetized 
Ba ion and the neutral Ba atoms streaming past it at time 100 m sec after 
burst and the differential energy distribution. There is much more energy 
in the peak near 4.2 km/sec than near 1-2 km/sec. but the collision frequency 
is much lower. Thus the more dense low velocity barium gas can couple 
energy from the neutral to the ions and magnetic field by collision processes 
while the high velocity gas cannot. The energy in the barium neutral gas 
travelling between 0 and 3 km/sec is 166 kJ, sufficient to create a magnetic 
bubble 680 m in diameter if the kinetic energy were totally expended. 

Obviously only a fraction of the energy of the neutral gas was transferred 
to the ions and magnetic field. Collisions are significant out to a few 
hundred meters g1 ven ^e assumpti on that^n^ly 1 5%_. of the bari urn was_ vapori zed . 

We have no firm data to verify the vapbration efficiency. If a greater 
percentage of the barium is vaporized the effect of neutral-ion coupling 
would be effective out to a larger radius. . 

The gases produced by the 4.5 kg of ^high explosives may .contribute 
to the snow plow effect sweeping out ions. The detonation of 4.5 kg of 
Comp B, would produce 141 moles of gas products, primarily N2» CO, CO2, CH3OH 
and lesser amounts of CH2O2, CkLj, NH3 and H2O. The gases at detonation have ■ 
a temperature of 4420*K and an expansion velocity of 7.84 km/sec (Cook, 1971). 
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If we assume all the gas Is In a spherical shell from 7 to 8 km/second, then 
the collision frequency between an atmospheric Ion, or a Ba Ion and the ex- 
plosion gas at time 100 m sec Is 1243/sec, and falls off as t"^. The 
actual density distribution's more likely to approach a Maxwellian distri- 
bution, and the collision frequency would be lower. Also the molecular 
weight of the gas products Is much less than 6a so the momentum transfer 
from the light gas molecules to Ba Ions will be less. The explosive gas 
products thus do not seem to be effective In sweeping out Ba Ions for more 
than the order of 100 m radius, and the velocity exceeds most of the barium 
gas, so most of the Ions would form behind the gas front. 

We also Investigated the possibility that the explosive gases might 
absorb solar UV, preventing Ionization of the expanding Ba gas for a period 
of time approximately 2.5 seconds. However, none of the gases has a 
significant absorption cross section In the 3000 A region, and the number 
of molecules In the path length Is at least four orders of magnitude below 
the amount necessary to produce significant UV absorption. 

Similarly the argument that the barium gas Is optically thick for 
about 2.5 seconds does not hold up to simple calculations. Also If the 
optical thickness were responsible for the black hole, the Inner edge would 
not have a steep gradient as the barium thfnned dut.~ Vet th^^ 
of the hole Is sharp (see Figure 4). . 

Barium Plasma Effects : Bubble Machine I 

The payload consisted of the 4.5 kg radial shaped charge on a spring loaded 
separable daughter section and an Instrumented section attached to the 
Tomahawk motor. All power to the Instrumented section and the separation 

t 

mechanism failed at lift off so only ground based optical and other obser- 
vations are available. Detonation occurred at 08 10.55 UT 25 March 1981 
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and at 66.08*’ lat. 146.47*’W long., 495.7 km altitude. Kp was 5- during 
the three hour interval of the experiment, but It was quiet magnetically 
at College during the event. As evidenced by the rapid convection of the 
Ba'*’ ions there was a substantial DC electric field. The event occurred 
well polewards of diffuse bands of aurora. The aurora was quiet and homo- 
geneous over the northern half of the sky. Meridian scanning photometer 
records In 6300 A, 5577 A and 4278 A indicate a characteristic energy 
deposition of approximately 600 eV to 1 KeV with a number flux of approxi- 
mately 10® electrons/cm^sec. The 4278 A emission was about 200 Rayleighs. 

Two photometers at Ester Dome, Alaska, with narrow bandwidth filters at 
4278 A were aimed near but below the burst point, and In the E region 
magnetically conjugate to the burst point to look for Induced auroral 
emissions. Figure 6 shows the traces of the two 4278 A photometers during 
the event. A damped oscillation of about 5 second period is clear from 
both records but there Is a delay of about 1 second at the 100 km level. 

The 500 km photometer was moved 3® in azimuth between 5 and 8 seconds 
after the event which had very little effect because of the high elevation 
angle. It was moved In elevation towards the neutral Ba cloud at 14 seconds, 
so the two traces are not expected to correlate after that time. The 
amplitude of the oscillation in 4278 A was somewhat larger in the 100 km 
trace. The auroral perturbation at both elevations begins with a decrease 
In the 4278 A background emission, and the modulation Is about 50%. The 
period of the oscillation increases with time: the first cycle Is about a 
4 second period, the second 5 seconds, the third about 6 seconds. 

Barium Plasma Effects : Bubble Machine II Bubble machine II was launched on 

March 15, 1982 with detonation at 06:43; 59.233 UT and at 66.13®N lat., 
146.325®W long, and 460.11 km altitude. . The Injection was viewed by use 
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Figure 6. Comparison of 5* field of view photometers with 4278A Interference 
filters looking near the Oubble Machine 1 burst region (upper 
trace) and In the 110 km region near the foot of the burst point 
field line. Mote the one second phase shift between the upper 
and lower traces. A damped oscillation In the 4270A emission 
with about a 5 second periodicity commences with a decrease. 

AURUS TOMAHAWK 34.005UE, BUBBLE MACHINE RADIAL SHAPED. CHARGE 

I 




of image orthicon TV systems at Ester Dome (64.88® N, 148.05® W) and Fort 
Yukon, Alaska (66.56® N, 145.22® W) and from the NASA Ames Lear jet flying 
over Prince George, British Columbia (53.81® N, 124.29® W) . The Learjet 
was positioned to provide as near as possible a side view of the injection. 
Figure 3 shows a view of the slow material (1-2 km/sec.) from the Learjet. 
The slow barium appears to have the form of an expanding doughnut. Figure 7 
shows a view of the neutral Ba gas at 5 seconds after detonation from Ft. 
Yukon looking almost perpendicular to the plane of the injection. The 
half angle of the doughnut out of the plane of injection measures 13®. 

This applies only to the slowest Ba gas as the TV was overdriven at earlier 
times when the fast jet material could have been seen. For calculations 
of Ba density and collision frequency we have assumed 13® for all the 
material. 

Although the magnetic activity at the time of detonation was quiet, a 
substantial E field existed in the detonation region. Although the classic 
King Crab black hole and stellate Ba ion structure was probably formed, the 
E X B drift rapidly distorted and blurred the structures. The observed 
average ion drift velocity was 1.06 km per second, implying average E 
fields: En = 38.9 mV/m and Ee * 2.8 mV/m. 

The main rocket payload with its Instrumentation could be observed on 
film and TV pictures. ‘ Triangulation and least mean square fitting to a 
trajectory has produced an estimate that at burst the mother was 0.384 km 
south, 0.505 west and 1.065 km below the burst point for a total s^eparati on 
of K24 km. The magnetic “field line threading through the instrumented 
payload was 824 m radially from the burst point in -the plane of the 
injection. The triangulated toroid dips about 2® to the northwest. 
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If we assume that the barium Ions were pushed out of an area within 
a 2.5 km radius as observed In the King Crab experiment (Wescott et a1., 

1980) In a period of 2.5 seconds as Implied from the slowest material, 
then we can build up a scenario as shown In Figure 8. The figure shows 
the location of the ring of Ions surrounding the black hole at +1, +2 and 
+3 seconds under the Influence of the background E field and the neutral 
collisions. The neutrals and ions would have originally passed over the 
Instrumented rocket In about 300 msec while the ionization process was 
in progress. Then for about 2 seconds the Instrumented rocket should have 
been within the Ba Ion depleted region. At about +2.5 seconds the E x B 
drift would have brought the concentrated ion ring back over the Instrument. 

The telluric E field at Poker Flat was measured using a 1.59 km grounded 
dipole oriented nearly magnetic N-S (at N 25*E). The data storage system 
was capable of recording 190 seconds of data taken at 3 samples per second. 
Figure 9 shows the telluric E field data which contains the time of the 
Bubble Machine II detonation. The data In Figure 9 have been smoothed using 
a three point running average. The time of detonation Is denoted as t=0 on 
the data plot and occurs approximately 20 seconds Into the data record. 

Two features stand above the background jjoj^se In the telluric data, 
the first Is a two-cycle wavelet commencing about 2 seconds after the burst 
and the second, similar wavelet of smaller amplitude occurs approximately 140 
seconds after detonation. Both of these pulses have periods of approximately 
5.5 seconds. A somewhat longer period variation occurs near the 100 second 
mark. The peak to peak amplitude of these pulses In the first wavelet Is 
about 10 mV/km. These are the only pulses of this type observed during a 
few hours data collection both before and after the experiment. 
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Figure 8. 
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^1"? surrounding the black hole of 
Bubble Machine II at one second Intervals. The position is due 

and^thP^?yR®rfr?fri® Ik P«sh the ions outwards to about 2.5 km 
and the EXB drift to the northwest. As a result the ions would 
sweep across the instrumented payload for about 300 ms then 
return starting about 2.5 seconds after burst. 
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close to magnetic north south at Poker Flat. Alaska, during the Bubble 
Machine II event. The two cycle pulsation commencing two sec. 
after burst; appears to be a result of the explosive Ba perturbation 
The two cycle wave commencing at 140 sec. Is consistent with 
AlfvCn wave' propagation to the southern conjugate ionosphere 
and return. 



The telluric system had been used In other experiments to measure the 
effect of auroral magnetic disturbances on telluric electric potentials. 
During these previous experiments a system response in the frequency 
range of Interest was established such that E » k dB/dt. For E in volt/km 
dB/dt in nT/sec, k » 0.1. (dB/dt is measured for the horizontal magnetic 
component perpendicular to the line.) The measured system response of 
10 mV/km would imply a wave in the east west component of the magnetic 
field of about 0.1 nT amplitude. 

Corresponding pulses could not be found in the magnetometer records 
from the rocket payload which had threshold sensitivity of 50 nT. Pulses 
were also absent in the induction coil magnetometer on the ground at College 
Alaska. The induction coil had a threshold noise level of 4x1 0"^ nT in the 
5 second period band. If the pulses are assumed to be electromagnetic 
wave estimates of the E/B ratio indicates the magnetic component ground 
amplitude would be approximately 1/10 the noise level and undetectable. 

We turn now to the results of the instrumented portion of the Bubble 
Machine II payload. A frequency-time-intensity (sonogram) plot of the VLF 
waves (10-16,000 Hz) detected during a time period centered on the burst is 
presented in Fig. 10. The data are from the electric field instrument. 

Natural -VLF emissions were detected-after the burst with frequency and 

temporal characteristics typical of similar auroral experiments (Temerin 
and Kelley, 1980). They display a sharp low frequency cut-off at the lower 
hybrid frequency. If we .assume that the undisturbed plasma is entirely 
composed of Ions, the observed lower hybrid frequency at a flight time 
of 360 seconds implies an ambient plasma density of 5 ± .4 x 10^ cm-3. 

There is no reason to suspect that these emissions were caused by the 
explosion. 
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Figure 10. A VLF sonogram showing electric field fluctuations intensity 
function of time and frequency. 
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The dark line at the burst time shows that Intense broadband emissions 
accompanied the explosion, with the signal exceeding the Instrument noise 
level up to at least 12,000 Hz. Low frequency emislons (^120 Hz) are also 
evident In the sonogram for several seconds (see discussion below) and an 
Isolated burst of mid-band signal centered at 4.3 kHz started about two ^ 
seconds after the burst. Directly after the burst the VLF signal seems to 
be suppressed for about one second at frequencies above the ambient lower 
hybrid frequency. 

The low frequency waves were also detected by a fixed bias Langmuir 

« 

probe showing their electrostatic nature. A sonogram of the low frequency 
(10-500 Hz) signal Is presented In the top panel of Fig. 11. The burst 
time Is Indicated and again the Intense Initial broad band noise pulse is 
apparent. This Is followed by a one second burst of signal at frequencies 
below 200 Hz. A gap In the emission strength then occurred for about two 
seconds which was followed by an extended period [-6 seconds) of emissions 
again well below 200 Hz. In fact, these waves display an extremely sharp 
high frequency cutoff just above 100 Hz. 

The second panel shows a signal proportional to logio of the current 
to a cylindrical fixed bias Langmuir probe. Care must be taken in Interpreting 
this output due to possible vehicle potential variations In the chemically 
active environment. We monitored the potential difference between the skin 
and one of the spheres. In fact, during the Initial post burst one second 
Interval this channel displayed a total peak to peak variation of about 0.4V 
with a waveform similar to the relative density waveform plotted in Fig. 11. 

We determined the value of dl/dV for the Langmuir probe during much of the 
flight (LaBelle et al., 1982) and hence can estimate the' apparent density 


representative power spectrum. 


March 15, 1982 ^ 




Poker Flat Alaska 







change under the assumption that the observed signal Is entirely due to the 
skin voltage change (which is the Langmuir probe reference potential). 

About one-half of the total response observed can be easily explained In 
this manner and we conclude that the Initial variations are not real density 
changes. However, by 0644:01 the vehicle potential returned to Its pre- 
explosion level and remained there for an extended period. We thus conclude 
that at 0644:02 the payload was Immersed In a plasma which was 25% more 
dense than the ambient plasma prior to the event. By 0644:11 the measured 
plasma density returned to the pre-explosion value. 

The second burst of low frequency waves Is well correlated with the 
density enhancement. They do not seem to be associated with the peak In the 
enhanced plasma density, but rather arrive on the “trailing edge" of the 
plasma. A similar relationship between emissions and plasma gradients has 
been observed In other chemical release experiments (Kintner et a1., 1980; 
Kintner and Kelley, 1982). Notice that the waves near 4 kHz also occurred 
at a time when the plasma density was decreasing. 

Representative («n/n)2 spectra are plotted In the bottom panel of 
Fig. 11 with labels corresponding to the time periods similarly numbered In 
the upper panel. The burst spectrum (#1) displays a (1/f)2 form. Fixed 
bias Langmuir probes are sensitive to electric field variations under some 
circumstances (Kelley et a1., 1982). However, comparison with the 6E 
spectrum shows that the «n/n variations are 15 db higher than those 
expected fran coupling to' the electric field. We thus conclude that 
electrostatic waves or a combination of electrostatic and electromagnetic 
waves are responsible for the burst waves. Similar conditions were met In 
considering burst related emissions In other experiments of this type. 



Spectrum 3 shows how quiet the region of peak plasma density was. For 
example, at 100 Hz the signal Is at least 30 db lower In spectrum 3 than In 
any of the others plotted. 

Spectrum 4 shows the extremely sharp nature of high frequency cut-off 
In the signal. The steep signal decay begins at 120 Hz and falls off 
roughly as f"^2 for f>120 Hz. Spectrum 2 also shows a high frequency cut- 
off above 10 Hz, but one which Is much more gradual. 

The Integrated fluctuation strength In the band 70 - 120 Hz Is plotted 
In Fig. 12 for the (6n/n) and (6E) receivers. The signal to noise 
ratio was better In the former detector due to the Interference bands in 
the electric field channel which can be seen at low frequencies In Fig. 10. 
The modulations In the signal are mirrored In the two channels and are not 
related In any obvious way to the spin rate. 

We are turning now to the dc electric field signature of the event. 
Previous explosive chemical release events have generated significant quasi- 
dc electric field probes. For example, the cesium lined spherical high 
explosive used In the Trigger experiment created a 200 mV/m electric field 
pulse (Holmgren et al., 1980). That pulse was followed several hundred ms 
later by an Intensive field aligned electron precipitation event, consistent 
with Instability of an upward propagating Alfven wave (Kelley et al . , 1980). 

In the present experiment, a quasl-dc electric field pulse was also observed 
but which lasted less than 100 ms and which had an amplitude of about 10 mV/m. 
A single dipole antenna can only measure one component of such a pulse and 
hence the measured value represents a lower limit on the field strength.. 
Nonetheless, It seems clear that the radial shaped charge was less efficient 
In generating an Alfven wave than the lower velocity (but higher plasma 
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density) cesium experiment. The magnetic field measurement on board the 

rocket was sensitive to perturbations greater than about 50 nT. For an 

Alfven speed of 500 km/s this corresponds to a 25 mV/m electric 

field perturbation. Thus, the absence of an observable magnetic field 

signal Is consistent with the small electric field pulse. (In the most 

recent more efficient radial shaped charge experiment conducted in April. 

1984 a much larger electric field perturbation was observed.) 

Figure 13 shows the Ion energy that was analyzed and the corresponding 

Ion differential energy flux as a function of time during nearly one second 

* 

around the burst time. 

The Ion ESA sweeps from 55eV to 10 keV In 125 msec, taking 12 steps of 

4 accumulation times each. The downsweep followed the same pattern, with 

the energy levels staggered to lay between the upsweep levels. The pitch 

0 

angle was within 10 degrees of the field. 

We group these data Into three periods of energetic precipitation: 
one centered at 6:43:59.33 UT. one long burst lasting between 6:43:59.40 
and 59.54. and one final Interval at 59.60. The energies observed ranged 
from the lowest level (nominally 55eV) to 6.8 keV (the highest energy level 
that the sweep Included was 10 keV). Since the highest energy particles pro- 
, duced In the explosion are -60 eV. kev energy particles must be assumed to 
be already present. The weak electric field (~10 mV/m) observed makes It 
unlikely that a large electrostatic shock could be responsible for accelerating 
the particles as was believed to be the case for Project Trigger (Holmgren. 

. 1980). ^ 

Although the time resoltulon Is not sufficient for a detailed dis- 
persion analysis. It Is worth noting that the observations do not fit a' 
classic dispersion pattern of arrival from a distant source. Indeed the first 
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particles that are observed at 59.33 are of comparatively lower energy. 

This makes It likely that the source of the particles was close to the payload. 
It may be that the Ions were already In existence and were simply scattered by 
ordinary Ion-Barium neutral collisions Into pitch angles near zero. This 
explanation would require that the distribution be sharply peaked at large 
pitch angles, since otherwise our detector (which was nominally aligned with 
the magnetic field) would have observed the energetic Ions prior to the release. 

Each period of Ion precipitation can then be correlated with the 
passage of each of the neutral rings (shown In Figure 8) through the field 
line connected with the Instrumental payload. 

The Ba expands radially In an Injection plane nearly perpendicular to 
the magnetic field. The time of arrival of the first ring at the field 
line connecting with the payload Is 100 msec after the explosion; the first 
energetic Ions are observed at 93 ± 18 msec after the explosion. This 
agreement Is excellent, although the transit time for a 0+ Ion of the 
appropirate energy adds another 10 msec or so to the 100 msec figure. It 
Is worth noting that the first Ions are observed prior to the direct time 
of flight of the fastest debris from the explosion to the detector. 

The second group of Ions Is first detected 172 ± 18 msec after the 
explosion, this corresponds very well” to'the time of arrival bf the faster 
atoms In the second Ba ring to the field line connecting with the payload. 

Such atoms are moving at ~4.4 km/sec (refer to figure 2). 

From figure 2 alone ^one might expect that the precipitation of Ions 
would be continuous (perhaps with Intensity variations) until after the 
passage of the slowest ring (peaked at 2 km/sec). Ions are Indeed observed 
until the sweep returns to lower energy levels (this Is why we Include the 
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entire time 6:43:59.405 to 59.546 as one burst). The end of this burst at 
59.546 corresponds to at transit velocity of 2.5 km/sec; there Is no 
particular matching feature on the velocity distribution chart. We speculate 
that this simply represents a depletion of available energetic Ions to be 
scattered; the passage of the previous 6a atoms could have by then essentially 
removed most of the energetic Ions. 

The final major observation of energetic Ions occurs at about 59.625; 
this coincides both with the arrival time of the 2 km/sec peak and with the 
detector reaching the lowest energy resolution level (nominally 55 ev). Of 
the two possible factors, the latter Is likely the most Important. This Is 
the only time near the explosion that the 55 ev level Is measured. It Is 
possible that there are far more comparatively low energy Ions available to 
be scattered that there are In the hundred ev range or higher. Indeed, the 
explosion Itself energized Ba atoms up to about 65 ev; some of these may 
have become Ionized or transferred energy to thermal plasma Ions. 

Figure 5 shows the estimated collision frequency of Ions with the 
neutral particles In the released cloud; It Is high enough (on the order of 
70/sec) that the scattering mechanism suggested here Is plausible. Barium 
(atomic weight J 37) Is of course much heavier th an the background plasma of _ 
mostly atomic ozygen and hydrogen, hence large angle scattering Is more 
likely than would otherwise be the case. 

The retarding potential analyzer, which was devoted to 240 eV electrons, 
recorded a single burst about 33 msec after the first Ions arrived, and 
lasting approximately 10 msec. These data are shown In Figure 14 along 
with the pitch angle of the analyzed particles. It Is not clear why there 
Is only this single electron event. The fixed energy was chosen to be at 
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the peak of the electron spectrum seen In the Trigger experiments (Lundin, 
1980). Thus, the absence of extended precipitation In these experiments 
contrasts greatly with those of Trigger. Both the particle data and the 
dc-electric results of Bubble Machine II imply a much weaker effect than In 
that experiment. 

In summary, we have Interpreted the particle results In terms of an 
already existing population of energetic Ions pitch angle scattered by the 
released barium Into the observation cone of the detector. This distribution 
must have been exclusively at large pitch angles In order to have avoided 
earlier detection by the ESA. Such a population Is called an Ion conic, 
and have been observed previously, although at higher atltltudes or lower 
energies In both (Gerney et al., 1981) and sounding rockets. 

Yau et al . (1983) observed transversely accelerated Ions (TAI) with pitch 
angles 90 to 130 In two sounding rocket flights into an active portion of 
the auroral zones at altitudes 400-600 km. These were In the range up to 
several hundred eV. The strongest TAI events tended to occur near regions 
of minimum energetic electron flux. 

Thus, at least In the case of Bubble Machine II, there Is little 
evidence that the energetic Ions were created by the burst. 
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Fiqure 13. The bottom panel shows the flux of Ions at the energy per charge 
shown In the bottom panel for about one second after the chemical 
Ignition. All of these Ions are at nearly 0 degree pitch angle. 
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Figure 14. The flux of electrons at 240 ev at the given pitch angle in the 
top panel . 



III. Discussion 

A. Stimulated and Auroral Pulsations and Earth Current Modulations 

The normal auroral luminosity vs height distribution has a maximum In 
the E layer. The altitude of the maximum depends upon the energy of the 
precipitating auroral electrons and the density of the upper atmosphere. 

The 4278A emission comes from the N 2 + first negative band. If the 
atmosphere above 400 km Is sunlit, resonance scattering from N 2 '*’ Ions 
Increases the volume emission rate so that the observed 4278A emission 
from above the sunlight line becomes significant with respect to the E 
layer peak. Figure 15 shows the estimated 4278A volume emission rate 
height profiles for Bubble Machine I. It Is clear that there was sufficient 
4278A emission near the burst where the photometer was pointed to allow 
the 501 modulation observed. Similarly at the lower photometer there is 
sufficient light to support a modulation. Deehr and Romick (1977) have 
reported a similar modulation of auroral emission due to a Ba thermal 
release near the same altitude. The question Is how did the plasma pertur- 
bation physically Induce the modulation In the light? 

Figure 16 shows the fields of view of the two photometers with respect 
to the detonation, the 30 km screening height, the geometrical umbra and 
the aurora. It Is conceivable that the two correlated auroral modulations 
are a coincidental natural pulsation In the main auroral luminosity near 
110 km, but the auroral wave would have to travel polewards at 210 km/second 
horizontally. Cresswell (1968) has reported fast auroral waves In pulsating 

I 

aurora with velocities ranging between 30 and 150 km/sec. The all sky TV 
located at Poker Flat does not show any detectable wave train, but the 
system Is not very sensitive. It seems more likely that an oscillating 
perturbation began near the burst point and propagated down the field line 

In the E-layer. This Implies a propagation velocity of 380 km/sec. 
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Figure 15. Volume emission rate of 4Z78A vs altitude. Bubble Machine I. 

^ Solid line Is the normal contribution of the precipitating 

electrons. The dashed curve Is the fluorescence scattering In 
the sunlit position. 

I 

I 

i 










OF POOR QUAdrrV 


Both the upper and lower photometer traces Indicate an initial intensity 
reduction followed by an increase over the background level. If the 
Injection initiated a magnetosonic wave train of a period of 5 sec., the 
first arrival at any point would be a compression of the magnetic field. An 
increase in the magnetic field would have the effect of raising the mirror 
height of the energetic electrons and vice versa, a following rarefaction 
of the field would lower the mirror height. This would also modulate the 
number of electrons precipitating and the volume emission rate. 

Using the theory of Lutomirski (1978) we have calculated that the wave 
Bubble Machine I produced should have an amplitude near the burst of 40 nl. 
This is the order of magnitude registered by the electric field detector 
during the Bubble Machine II experiment of the following year. Assuming a 
50 nT wave amplitude, the mirror height could be raised or lowered by 2.4 km, 
or an overall change of 4.8 km at the E layer. 

Thus if modulation of the 4278A light is due to a magnetosonic wave, 
a larger amplitude wave than was measured in situ is suggested. Since the 
electric and magnetic field detectors were below the level of the burst 
plane and about 824 m from the field line intersection, it may be that the 

-geometry reduced the amplitude of the aB at the mother payload 

Another possibility is that the electric field component of the 5 
second period wavetrain has a small component parallel to B which might be 
effective in modulating the precipitating electrons. 

J ■ ■■ ' 

B. In Situ Wave Observations on Bubble Machine II ' 

We have already noted that the quasi-dc electric field signal implies 
that a modest amplitude Alfven wave was produced in the expansion. 

As discussed below, the amplitude is not inconsistent with estimates based 
on the ground wave measurements. 
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The higher frequency waves shown In Figs. 10-12 show a good correlation 
with the optical data. For example, the signals below 100 Hz appear in two 
sets— one lasting about one second and the other abut 8 seconds— with 
about a two second gap between the emissions. Referring to Fig. 8, we can 
understand the emissions sequences as follows. The burst waves are associated 
with the explosion and leading edge of the Ion ring. The low frequency 
waves 120 Hz) probably lie on the density gradient. The gap in signal 
Is totally consistent with the black hole since it should be absent of any 
barium related emissions. When the ion ring convects back across the 
payload some 2.5 seconds later, we again see the low frequency waves, this 
time clearly associated with the decreasing density gradient. The ratio of 
the durations of the two low frequency wave events are also consistent with 
an expansion/convection model associated with the ion ring. 

As discussed in the next section, the burst waves may play a role in 
the black hole formation process. Here we investigate the lower frequency 
wave emission mechanism (f<120 Hz). The barium plasma feels a radially 
outward force of the form 

F = mvi^ 

where £n Is the neutral _expans ion . velocijty and.vin is the ion neutral 

collision frequency. This force is equivalent to a gravitational acceleration 
given by £' ■ u.n which is subject to the Rayleigh-Taylor (R-T) 
instability in regions of inward plasma density gradient (£' must be anti- 
parallel to £n). ' ' ■ 

-The growth fate for the collisionless R-T process is 

y o (g7L)l/2 ^ i\) 

while for the colli si onal R-T process it is 

y n g'/vinL (2) 
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where L Is the gradient scale length. The bounding relationship 
between the colllslonless versus the colHslonal R-T process (Ott, 1978) 

Is v^n* ■ 4u/L » 160 Hz which, according to Fig. 5 Is greater than v-|n 
at 100 BIS. A transition thus occurs from a colllslonal process at early 
time to a colllslonless process later. At t>1 sec, the optical data Indicate 
L • 100 m while at .1 seconds It should be considerably smaller. If we let 
L»40 m and u,, ■ 4000 m/s and take v,jp * 75 sec"^ (from fig. 5), we find 
Y > 80 Hz which implies an e folding time of about 12 ms. Since the neutral 
gas Is observed to be structured, there will certainly be Initial finite 
perturbations to feed upon this energy source. The waves associated with 
this process should peak at kL • 1 which Implies X - 300 m. If we Interpret 
the frequencies In Fig. 11 (spectrum 2) as Doppler shifted zero frequency 
turbulence (R-T waves have very low phase velocity), the break in the spectrum 
to a power law occurs at about 50 Hz which corresponds to wavelengths greater 
than about 100 m If the waves are being carried out at the expansion speed. 
Fluo-tuatlons are detected to about 1200 Hz which Is roughly 3 m.( under this 
same zero frequency hypothesis). 

The wave spectrum number 4 Is also quite Interesting. First, the 
waves are located on the negative density gradient which corresponds to the 
unstable "leading edge" In the expansion phase. ^Second, “this edge 1s“ “ ~ 

further unstable to the ^ x ^ Instability (^ x _B parallel to vn) as the 
structure convects back toward the northwest. At 1.06 km/s velocity, the 
very steep decrease In signal at 120 Hz corresponds to about 8 m wavelength. 
This Is the gyro radius of a 1400" K barium Ion In^a .5 Gauss magnetic field. 

The cutoff Is so sharp that some characteristic frequency or scale 
must be Involved. The Ion gyro frequencies are all well below 120 Hz and 
the lower hybrid frequencies are all well above. It thus seems likely 


that the 8 m barium gyro radius (r^) is the important parameter. Furthermore, 
it is wen known (Coppi, 1965) that finite ion gyro radius (FLR) effects stabili 
the generalized Rayleigh-Taylor instability. A case can thus be made that 
the observed plasma structuring is due to a flute-mode Rayleigh-Taylor 
process initiated by the neutral gas expansion and perhaps sustained by 
the £ X Bi instability after the expansion ceases. 

Pursuing this line of argument further we may have tested a recent 
theoretical prediction by Sperling and Krall (1981) who investigated the 
"prompt striations" observed in a shaped charged experiment by Simons et al., 
(1980). They showed that the FLR stabilization only occurs if the ratio of 
background plasma density to the barium plasma density is above a certain 
threshold value. They argued then that in a low density background, short 
wavelength waves (X • r<) should occur early in the expansion process, 
creating the observed prompt striations. Our data is in excellent agreement 
with this proposal. For example, in the initial phase (spectrum 2) waves 
are present with a power law form between 3 m and a few hundred meters. Note 
that the gyro-radius associated with the beam velocity rather than the thermal 
velocity, is about 40 m. On the other hand when the expanding ring convects 
back to the payload we find no such power law and in fact an extremely sharp 
cut-off near (the thermal value). 

From Figures la-c of Sperling and Krall (1981) for the returning plasma, 
which has their parameter 7/Pi * 0.1, we can conclude quantitatively 
the short wavelength wavds should be totally stabilized as observed. For the 
early phase it is more difficult to determine p/pi, the ratio of ambient 
mass density to beam mass density, but it certainly can only be a lower value 
which is consistent with the theory and the observed existence of short wave- 
length waves. 
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C. Black Hole Formation 

The burst waves seem to be a mixture of electromagnetic and electrostatic 
waves with frequencies up to and in excess of the lower hybrid frequency. 

These may play a role in the black hole formation If the waves provide an 
anomalous coupling between the barium ions and the neutral gas. This can 
be seen as follows. The ion equation of motion is 

nM^dV^/dt « neE ♦ ne(^f x B) + nM^v^n Up 
If the ions will follow the neutral gas and be swept out of 

the region near the explosion point. Now, classical collisions fall 
below ni within a few hundred ms. However, waves may destroy any organized 
gyro motion and allow efficient coupling to the neutrals. The waves 
themselves may be driven unstable by the modified two stream instability 
caused by the £ x £ streaming of electrons subject to the electric field 
pulse from the explosion. The threshold for this process would be about 
40 mV/m for the conditions of the experiment. 

The polarization E field due to the large gyroradius of the Ba+ ion 
with respect to the electrons and the coupling of this field to the conducting 
ionosphere below must be considered. The polarization field in the first 

approximation is radially inwards. .It would propagate along the magnetic 

field at the Alfven velocity. The Alfven is estimated to be about 
2000 km/sec at a nighttime electron density of 10^ electrons/cc. Thus 
the polarization field could be shorted out in the ionosphere with a time 

I 

4elay of the order of 200 m sec. The field should have caused an upwards 
current near the burst point field line and downwards cylindrical current 
near the outer edge of the barium gas. 

The current should primarily be carried by the electrons, but the net 
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effect would also be to accelerate Ba+ away front the plane of the explosion 
In a cylinder expanding In radius with time. The TV observations cannot 
confirm definitely the parallel to B motion due to the overwhelming bright- 
ness of the Ba neutral gas in the line of sight to the field line. 

In Bubble Machine II the Cornell AC E field dipole measured some CH 
lower hybrid resonance electrostatic waves up to about 500 ni sec subsequently 
followed by a few seconds of waves at harmonics of the 0+ gyrof requency . 

The upper frequency limit of the AC detector was 15 kHz, so It was not 
possible to observe the occurrence of upper hybrid waves. The suggestion 
Is made that small scale turbulence can be effective In allowing plasma 
transport across magnetic field lines. 

The total effects of the momentum coupling of the neutral Ba gas and 
the explosive gas products to the Ba**- Ions, the polarization field and 
coupling to the Ionosphere with a time delay, may allow the ions to move 
outward freely to the 2.5 km radius required, or electrostatic turbulence 
may allow the ion motion. Self-consistent calculations on a big computer 
are required to answer the questions of the cause of black hole. 

D. Discussion of E-Fleld Pulsations 

It is tempting to identify the second wavelet as an echo of the original 
burst wavelet which has travelled along the field line to the opposite 
hemisphere and back. Periodic emissions In the Pci (0.3-3Hz) band are 
coninon at College, Alaska and exhibit repitlon times from two to five 
minutes. These are Alfveh waves which travel through the plasma In 
the ion cyclotron mode with a group velocity which depends upon the ion 
density of the magnetospherlc plasma. The 140 second delay and the 5.5 
second period observed In the telluric record are compatible with wave 
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patterns c<»nmonly observed on the L*5.5 shell at College. However, since 
no magnetic signature was observed, either at the rocket or at the ground, 
we cannot be certain that the detected pulses were electromagnetic. It is 
possible they represent electrostatic responses of the F-reglon wave guide 
to the burst. 

It Is, of course, possible that the wavelets are naturally occuring 
fluctuations detected by the telluric system. Since the recording system 
could only store 190 seconds of data we cannot check the data to see if 
similar pulses occurred 136 seconds before or 278 seconds afterwards. The 
wavelets at 2 and 140 seconds are the largest variations In the recording 
during the recording Interval. The fact that they occur after the burst 
and that they are compatible with the common period and propagation ranges 
of waves observed at College give support to the Interpretation that they 
were produced by the burst. Further, the observation of an oscillation of 
essentially the same five second period In the 4278A auroral emission 
following Bubble Machine I strengthens the case that these wavelets are 
burst-induced. 

The burst occurred at an altitude of 460 km. One minute after the 
burst an lonogram was taken which showed spread-F features wHh foF of 4.2 
MHz at 350 km with foFE_pf 3.5_MH2_at IIP km and no P layer. I f we assume 
simple Alfven wave propagation between the burst point and the E layer 
the lonogram allows us to estimate an Alfvdn velocity ~ 230 km/sec 
and hence a propagation time of about 1.4 seconds from the burst point to 
'the E region, this Is less than the observed two second delay between the 
burst and the first wavelet In Figure 9. However It Is of the same order 
as the time delay observed between the burst and E region. response observed 
on Bubble Machine I. 
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Another effect to be considered Is the size of the Ion depleted region 
and the velocity of the slowest Ba gas. The black hole radius Is about 
2.5 km and the lowest gas velocity about 1 km/second. Thus It Is plausible 
that the kinetic energy of the Ba gas Is effective In compressing the 
magnetic field controlled plasma for about 2.5 seconds. The relaxation of 
the plasma and overshoot Is consistent with a 5 second period oscillation 
continuing for several cycles before being damped out. However, the theory 
of Lutomirski (1968) would Indicate that the slow moving plasma would be very 
inefficient In radiating hydromagnetic waves. 

IV. Summary and Conclusions 

The radial shaped charge experiments discussed In this paper made 
waves and created some other curious effects. Due to circumstances and 
limitations in the range of equipment fielded for each experiment effects 
were not observed on all three experiments. The formation of the black 
hole was best observed on King Crab because of the absence of any convection 
electric field to distort the Ion configuration. However, stimulated 
auroral emissions and waves were masked by natural background noise at the 
burst time. The long period (5 second) modulation of 4278A auroral 
emissions was only observed following Bubble Machine I, because no photometers 

were deployed during Bubble Machine ii^. — Conversely the E field wave of 

period 5.5 seconds In telluric currents was only observed following Bubble 
Machine II, because no recordings were attempted during the other two 
experiments. ' '' ‘ 

Examination of possible physical mechanisms for the creation of the 
Ion-depleted black hole of radius about 2.5 km has not produced a viable 
explanation. Momentum transfer from the neutral gas to the Ions can probably 
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move the Ions out to a radius of about 500 m. Some other process allowed 
the Ions to walk out another 2 km as If there were no magnetic field. 

There Is insufficient energy In the gas to create a magnetic bubble that 
large, and the magnetic signatures Is at most 50 nT In situ. 

A possibility yet to be tested In a computer simulation Is that the 

t. : 

electron Hall drift was sufficient to set up a modified two stream Instability 
producing small-scale turbulence. The high frequency waves might allow the 
plasma to walk out the additional two kms. 

The Injection of the plasma radially apparently excited an eigen- 
frequency wave of period about 5 sec. In Bubble Machine I this was evidenced 
by modulation of the electron precipitatin causing 4278A auroral emissions. 
Based on a one second delay between the modulation near the burst, and In 
the E region at the foot of the field line, the propagation velocity was 
about 380 km/sec, within the Alfvdn propagation range. The magnetic 
pertubatlon does not seem large enough to raise or lower the mirror height 
sufficiently to account for the 501 modulation In luminosity. We suggest 
that a small E field component of the wave, parallel to Be, could more 
efficiently modulate the flow of electrons. 

The second manifestation of the 5 second period wave was In telluric 
field recordings during Bubble Machine II. The ^o cycle wavelet conmenced 
2 seconds after the burst, and was followed 138 seconds later by what we 
conclude Is probably a wave reflected frcxn the southern hemisphere Ionosphere 

I 

with about 1/2 amplitude. The Ionospheric density deduced from an lonogram 
taken one minute after detonation would predict a burst to E layer propagation 
of 230 km/sec, which would account for a time delay of 1.4 seconds Instead 
of the two seconds observed. The frequency. Interhemisphere bounce time 

46 













and F to E layer propagation velocity are well within observed values of 
natural waves at L ■ 5.5 for 1on«cyclotron mode waves. 
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FIGURE CAPTIONS 
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Figure 3. 
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Figure 6. 
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Figure 8. 


Cross section schematic drawing of radial high explosive shaped 
charge used in Bubble Machine I and II. Detonation is at the 
center causing the liner to collapse and Ba gas to jet out 
horizontally with radial syimetry. 

Differential 6a atom numbers vs radial velocity distribution. 
Bubble Machine II, derived from TV signal level scans. Fifteen 
per cent vaporiztion of liner was assumed to assign number scale 
value. 

Side view of Bubble Machine II from geophysical Institute TV on 
NASA Ames Learjet, 21. S sec. after detonation. The plane of 
the low velocity Ba gas doughnut is 17* with respect to the 
viewer and the distance is 1950 km. The two stars are v Andromeda 
and B.S. 157. 

TV view of the stellate ion structure produced in the March 16, 
1980 King Crab experiment 205 seconds after detonation. The 
inner hub region depleted of ionization is called a black hole. 

An almost complete absence of an ionospheric E field allowed the 
structure to be seen for about 20 minutes. 

Plots of the differential kinetic energy vs velocity. Bubble 
Machine II experiment and the collision frequency between a barium 
ion and the streaming neutral barium gas at 100 msec after the 
burst. Note that although the highest energy is centered near 
4 km/sec., collisions are too infrequent to couple energy from 
the neutral gas to the ions and magnetic field. Energy can be 
transferred near the low velocity peak near 1-2 km/sec. 

Comparison of 5* field of view photometers with 4278A interference 
filters looking near the Bubble Machine I burst region (upper 
trace) and in the 110 km region near the foot of the burst point 
field line. Note the one second phase shift between the upper 
and lower traces. A damped oscillation in the 4278A emission 
with about a 5 second periodicity commences with a decrease. 


Geophysical Institute TV picture of Bubble Machine II at *•■5 seconds 
from Fort Yukon, Alaska. Note the regions of blob-like irregularities. 

Plan view of part of the ion ring surrounding the black hole of 
Bubble Machine II at one second intervals. The position is due 
to the forces .acting to push the ions outwards to about 2.5 km 
and the EXB drift to the northwest. As a result the ions would 
sweep across the Instrumented payload for about 300 ns then 
return starting about 2.5 seconds after burst. 
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Figure 9. Recording of the voltage across a 1.6 km grounded dipole oriented 

close to magnetic north south at Poker Flat, Alaska, during the Bubble 
Machine II event. The two cycle pulsation conmenclng two sec. 
after burst appears to be a result of the explosive Ba perturbation 
The two cycle wave conmenclng at 140 sec. Is consistent with 
Alfvdn wave propagation to the southern conjugate Ionosphere 
and return. 

Figure 10. A VLF sonogram showing electric field fluctuations Intensity as a 
function of time and frequency. 

Figure 11. A low frequency sonogram showing the low frequency waves and 
representative power spectrum. 

Figure 12. Integrated electric field and density fluctuation signal strengths 
versus time. 

Figure 13. The bottom panel shows the flux of Ions at the energy per charge 
shown In the bottom panel for about one second after the chemical 
Ignition. All of these Ions. are at nearly 0 degree pitch angle. 

Figure 14. The flux of electrons at 240 ev at the given pitch angle In the 
top panel . 

Figure IS. Volume emission rate of 4278A vs altitude. Bubble Machine I. 

Solid line Is the normal contribution .of the precipitating 
electrons. The dashed curve Is the fluorescence scattering In 
the sunlit position. ^ 

Figure 16. Side view of magnetic meridian showing the location of the Bubble 
Machine 1 burst point, the magnetic field lines, the fields of 
view of the two 4278A photometers from Ester Dome, Alaska, and the 
aurora. 



